Gas-Phase Oxidation of Atmospherically Relevant Unsaturated Hydrocarbons by Acyl Peroxy Radicals.

This study assesses the atmospheric impact of reactions between unsaturated hydrocarbons such as isoprene and monoterpenes and peroxy radicals containing various functional groups. We find that reactions between alkenes and acyl peroxy radicals have reaction rates high enough to be feasible in the atmosphere and lead to high molar mass accretion products. Moreover, the reaction between unsaturated hydrocarbons and acyl peroxy radicals leads to an alkyl radical, to which molecular oxygen rapidly adds. This finding is confirmed by both theoretical calculations and experiments. The formed perester peroxy radical may either undergo further H-shift reactions or react bimolecularly. The multifunctional oxygenated compounds formed through acyl peroxy radical + alkene reactions are potentially important contributors to particle formation and growth. Thus, acyl peroxy radical-initiated oxidation chemistry may need to be included in atmospheric models.

Bicyclic guanidine superbase carboxylate salts for cellulose dissolution.

Bicyclic guanidines are utilized in organic synthesis as base catalysts or reagents. They also offer a platform for coordination chemistry, for example in CO2 activation, and their carboxylate salts offer an efficient media for cellulose dissolution. We have studied a series of bicyclic guanidines with varying ring sizes and with varying methyl substituents with a specific aim to find hydrolytically stable acetate salts for dissolution and processing of cellulose. Different superbase synthesis pathways were tested, followed by hydrolytic stability and cellulose dissolution capacity tests. The synthesis pathways were designed to enable the scale up of the production of the superbases considering the availability of the starting molecules and the feasibility of the synthesis. As a result, we found several hydrolytically stable bicyclic guanidine structures, which can overcome many of the reoccurring problems as carboxylate salts or free bases.

Multireference and Coupled-Cluster Study of Dimethyltetroxide (MeO4Me) Formation and Decomposition.

Peroxyl radicals (RO2) are important intermediates in the atmospheric oxidation processes. The RO2 can react with other RO2 to form reactive intermediates known as tetroxides, RO4R. The reaction mechanisms of RO4R formation and its various decomposition channels have been investigated in multiple computational studies, but previous approaches have not been able to provide a unified methodology that is able to connect all relevant reactions together. An apparent difficulty in modeling the RO4R formation and its decomposition is the involvement of open-shell singlet electronic states along the reaction pathway. Modeling such electronic states requires multireference (MR) methods, which we use in the present study to investigate in detail a model reaction of MeO2 + MeO2 → MeO4Me, and its decomposition, MeO4Me → MeO + O2 + MeO, as well as the two-body product complexes MeO···O2 + MeO and MeO···MeO + O2. The used MR methods are benchmarked against relative energies of MeO2 + MeO2, MeO4Me, and MeO + MeO + O2, calculated with CCSD(T)/CBS, W2X, and W3X-L methods. We found that the calculated relative energies of the overall MeO2 + MeO2 → MeO4Me → MeO + O2 + MeO reaction are very sensitive to the chosen MR method and that the CASPT2(22e,14o)-IPEA method is able to reproduce the relative energies obtained by the various coupled-cluster methods. Furthermore, CASPT2(22e,14o)-IPEA and W3X-L results show that the MeO···O2 product complex + MeO is lower in energy than the MeO···MeO complex + O2. The formation of MeO4Me is exothermic, and its decomposition is endothermic, relative to the tetroxide. Both the formation and decomposition reactions appear to follow pathways with no saddle points. According to potential energy surface scans of the decomposition pathway, the concerted cleavage of both MeO···O bonds in MeO4Me is energetically preferred over the corresponding sequential decomposition.

Effect of an Electric Field on the Structure and Stability of Atmospheric Clusters.

We study the influence of an applied electric field on the structure and stability of some common bimolecular clusters that are found in the atmosphere. These clusters play an important role in new particle formation (NPF). For low values of the electric field (i.e., |E| ≤ 0.01 V Å-1), we demonstrate that the field response of the clusters can be predicted from simply calculating the dipole moment of the cluster and the dipole moments of the constituent molecules and that the influence on the association energy of the cluster is minimal (i.e., <0.5 kcal mol-1). For higher field strengths |E| > 0.2 V Å-1, there can be more dramatic effects on both structure and energetics, as the induced dipole, charge transfer, and geometric distortion play a larger role. Although such large fields are not very relevant in the atmosphere, they do exist in some situations of experimental interest, such as near interfaces and in intense laser fields.

Computational Tools for Handling Molecular Clusters: Configurational Sampling, Storage, Analysis, and Machine Learning.

Computational modeling of atmospheric molecular clusters requires a comprehensive understanding of their complex configurational spaces, interaction patterns, stabilities against fragmentation, and even dynamic behaviors. To address these needs, we introduce the Jammy Key framework, a collection of automated scripts that facilitate and streamline molecular cluster modeling workflows. Jammy Key handles file manipulations between varieties of integrated third-party programs. The framework is divided into three main functionalities: (1) Jammy Key for configurational sampling (JKCS) to perform systematic configurational sampling of molecular clusters, (2) Jammy Key for quantum chemistry (JKQC) to analyze commonly used quantum chemistry output files and facilitate database construction, handling, and analysis, and (3) Jammy Key for machine learning (JKML) to manage machine learning methods in optimizing molecular cluster modeling. This automation and machine learning utilization significantly reduces manual labor, greatly speeds up the search for molecular cluster configurations, and thus increases the number of systems that can be studied. Following the example of the Atmospheric Cluster Database (ACDB) of Elm (ACS Omega, 4, 10965-10984, 2019), the molecular clusters modeled in our group using the Jammy Key framework have been stored in an improved online GitHub repository named ACDB 2.0. In this work, we present the Jammy Key package alongside its assorted applications, which underline its versatility. Using several illustrative examples, we discuss how to choose appropriate combinations of methodologies for treating particular cluster types, including reactive, multicomponent, charged, or radical clusters, as well as clusters containing flexible or multiconformer monomers or heavy atoms. Finally, we present a detailed example of using the tools for atmospheric acid-base clusters.

Correction: Ring-opening yields and auto-oxidation rates of the resulting peroxy radicals from OH-oxidation of α-pinene and ß-pinene.

[This corrects the article DOI: 10.1039/D2EA00133K.].

Molecular rearrangement of bicyclic peroxy radicals is a key route to aerosol from aromatics.

The oxidation of aromatics contributes significantly to the formation of atmospheric aerosol. Using toluene as an example, we demonstrate the existence of a molecular rearrangement channel in the oxidation mechanism. Based on both flow reactor experiments and quantum chemical calculations, we show that the bicyclic peroxy radicals (BPRs) formed in OH-initiated aromatic oxidation are much less stable than previously thought, and in the case of the toluene derived ipso-BPRs, lead to aerosol-forming low-volatility products with up to 9 oxygen atoms on sub-second timescales. Similar results are predicted for ipso-BPRs formed from many other aromatic compounds. This reaction class is likely a key route for atmospheric aerosol formation, and including the molecular rearrangement of BPRs may be vital for accurate chemical modeling of the atmosphere.

Influence of plasmons on the luminescence properties of solvatochromic merocyanine dyes with different solvatochromism.

The effect of localized surface plasmon resonance (LSPR) of a system consisting of a highly dipolar merocyanine dye and a silver nanoparticle (NP) was studied experimentally and theoretically. A theoretical model for estimating the fluorescence quantum yield (φfl) using quantum chemical calculations of intramolecular and intermolecular electronic transition rate constants was developed. Calculations show that the main deactivation channels of the lowest excited singlet state of the studied merocyanines are internal conversion (kIC(S1 → S0)) and fluorescence (kr(S1 → S0)). The intersystem-crossing transition has a low probability due to the large energy difference between the singlet and triplet levels. In the presence of plasmonic NPs, the fluorescence quantum yield is increased by a factor of two according to both experiment and computations. The calculated values of φfl, when considering changes in kr(S1 → S0) and the energy-transfer rate constant (ktransfer) from the dye to the NP was also twice as large at distances of 6-8 nm between the NP and the dye molecule. We also found that the LSPR effect can be increased or decreased depending on the value of the dielectric constant (εm) of the environment.

Atomic structures, conformers and thermodynamic properties of 32k atmospheric molecules.

Low-volatile organic compounds (LVOCs) drive key atmospheric processes, such as new particle formation (NPF) and growth. Machine learning tools can accelerate studies of these phenomena, but extensive and versatile LVOC datasets relevant for the atmospheric research community are lacking. We present the GeckoQ dataset with atomic structures of 31,637 atmospherically relevant molecules resulting from the oxidation of α-pinene, toluene and decane. For each molecule, we performed comprehensive conformer sampling with the COSMOconf program and calculated thermodynamic properties with density functional theory (DFT) using the Conductor-like Screening Model (COSMO). Our dataset contains the geometries of the 7 Mio. conformers we found and their corresponding structural and thermodynamic properties, including saturation vapor pressures (pSat), chemical potentials and free energies. The pSat were compared to values calculated with the group contribution method SIMPOL. To validate the dataset, we explored the relationship between structural and thermodynamic properties, and then demonstrated a first machine-learning application with Gaussian process regression.

Large Gas-Phase Source of Esters and Other Accretion Products in the Atmosphere.

Dimeric accretion products have been observed both in atmospheric aerosol particles and in the gas phase. With their low volatilities, they are key contributors to the formation of new aerosol particles, acting as seeds for more volatile organic vapors to partition onto. Many particle-phase accretion products have been identified as esters. Various gas- and particle-phase formation pathways have been suggested for them, yet evidence remains inconclusive. In contrast, peroxide accretion products have been shown to form via gas-phase peroxy radical (RO2) cross reactions. Here, we show that these reactions can also be a major source of esters and other types of accretion products. We studied α-pinene ozonolysis using state-of-the-art chemical ionization mass spectrometry together with different isotopic labeling approaches and quantum chemical calculations, finding strong evidence for fast radical isomerization before accretion. Specifically, this isomerization seems to happen within the intermediate complex of two alkoxy (RO) radicals, which generally determines the branching of all RO2-RO2 reactions. Accretion products are formed when the radicals in the complex recombine. We found that RO with suitable structures can undergo extremely rapid C-C ß scissions before recombination, often resulting in ester products. We also found evidence of this previously overlooked RO2-RO2 reaction pathway forming alkyl accretion products and speculate that some earlier peroxide identifications may in fact be hemiacetals or ethers. Our findings help answer several outstanding questions on the sources of accretion products in organic aerosol and bridge our knowledge of the gas phase formation and particle phase detection of accretion products. As esters are inherently more stable than peroxides, this also impacts their further reactivity in the aerosol.

Computational Investigation of Substituent Effects on the Alcohol + Carbonyl Channel of Peroxy Radical Self- and Cross-Reactions.

Organic peroxy radicals (RO2) are key intermediates in atmospheric chemistry and can undergo a large variety of both uni- and bimolecular reactions. One of the least understood reaction classes of RO2 are their self- and cross-reactions: RO2 + R'O2. In our previous work, we have investigated how RO2 + R'O2 reactions can lead to the formation of ROOR' accretion products through intersystem crossings and subsequent recombination of a triplet intermediate complex 3(RO···OR'). Accretion products can potentially have very low saturation vapor pressures, and may therefore participate in the formation of aerosol particles. In this work, we investigate the competing H-shift channel, which leads to the formation of more volatile carbonyl and alcohol products. This is one of the main, and sometimes the dominant, RO2 + R'O2 reaction channels for small RO2. We investigate how substituents (R and R' groups) affect the H-shift barriers and rates for a set of 3(RO···OR') complexes. The variation in barrier heights and rates is found to be surprisingly small, and most computed H-shift rates are fast: around 108-109 s-1. We find that the barrier height is affected by three competing factors: (1) the weakening of the breaking C-H bond due to interactions with adjacent functional groups; (2) the overall binding energy of the 3(RO···OR'), which tends to increase the barrier height; and (3) the thermodynamic stability of the reaction products. We also calculated intersystem crossing rate coefficients (ISC) for the same systems and found that most of them were of the same order of magnitude as the H-shift rates. This suggests that both studied channels are competitive for small and medium-sized RO2. However, for complex enough R or R' groups, the binding energy effect may render the H-shift channel uncompetitive with intersystem crossings (and thus ROOR' formation), as the rate of the latter, while variable, seems to be largely independent of system size. This may help explain the experimental observation that accretion product formation becomes highly effective for large and multifunctional RO2.

An Experimental and Master Equation Investigation of Kinetics of the CH2OO + RCN Reactions (R = H, CH3, C2H5) and Their Atmospheric Relevance.

We have performed direct kinetic measurements of the CH2OO + RCN reactions (R = H, CH3, C2H5) in the temperature range 233-360 K and pressure range 10-250 Torr using time-resolved UV-absorption spectroscopy. We have utilized a new photolytic precursor, chloroiodomethane (CH2ICl), whose photolysis at 193 nm in the presence of O2 produces CH2OO. Observed bimolecular rate coefficients for CH2OO + HCN, CH2OO + CH3CN, and CH2OO + C2H5CN reactions at 296 K are (2.22 ± 0.65) × 10-14 cm3 molecule-1 s-1, (1.02 ± 0.10) × 10-14 cm3 molecule-1 s-1, and (2.55 ± 0.13) × 10-14 cm3 molecule-1 s-1, respectively, suggesting that reaction with CH2OO is a potential atmospheric degradation pathway for nitriles. All the reactions have negligible temperature and pressure dependence in the studied regions. Quantum chemical calculations (ωB97X-D/aug-cc-pVTZ optimization with CCSD(T)-F12a/VDZ-F12 electronic energy correction) of the CH2OO + RCN reactions indicate that the barrierless lowest-energy reaction path leads to a ring closure, resulting in the formation of a 1,2,4-dioxazole compound. Master equation modeling results suggest that following the ring closure, chemical activation in the case of CH2OO + HCN and CH2OO + CH3CN reactions leads to a rapid decomposition of 1,2,4-dioxazole into a CH2O + RNCO pair, or by a rearrangement, into a formyl amide (RC(O)NHC(O)H), followed by decomposition into CO and an imidic acid (RC(NH)OH). The 1,2,4-dioxazole, the CH2O + RNCO pair, and the CO + RC(NH)OH pair are atmospherically significant end products to varying degrees.

Critical Role of Iodous Acid in Neutral Iodine Oxoacid Nucleation.

Nucleation of neutral iodine particles has recently been found to involve both iodic acid (HIO3) and iodous acid (HIO2). However, the precise role of HIO2 in iodine oxoacid nucleation remains unclear. Herein, we probe such a role by investigating the cluster formation mechanisms and kinetics of (HIO3)m(HIO2)n (m = 0-4, n = 0-4) clusters with quantum chemical calculations and atmospheric cluster dynamics modeling. When compared with HIO3, we find that HIO2 binds more strongly with HIO3 and also more strongly with HIO2. After accounting for ambient vapor concentrations, the fastest nucleation rate is predicted for mixed HIO3-HIO2 clusters rather than for pure HIO3 or HIO2 ones. Our calculations reveal that the strong binding results from HIO2 exhibiting a base behavior (accepting a proton from HIO3) and forming stronger halogen bonds. Moreover, the binding energies of (HIO3)m(HIO2)n clusters show a far more tolerant choice of growth paths when compared with the strict stoichiometry required for sulfuric acid-base nucleation. Our predicted cluster formation rates and dimer concentrations are acceptably consistent with those measured by the Cosmic Leaving Outdoor Droplets (CLOUD) experiment. This study suggests that HIO2 could facilitate the nucleation of other acids beyond HIO3 in regions where base vapors such as ammonia or amines are scarce.

Isomer-Resolved Mobility-Mass Analysis of α-Pinene Ozonolysis Products.

Highly oxygenated organic molecules (HOMs) are important sources of atmospheric aerosols. Resolving the molecular-level formation mechanisms of these HOMs from freshly emitted hydrocarbons improves the understanding of aerosol properties and their influence on the climate. In this study, we measure the electrical mobility and mass-to-charge ratio of α-pinene oxidation products using a secondary electrospray-differential mobility analyzer-mass spectrometer (SESI-DMA-MS). The mass-mobility spectrum of the oxidation products is measured with seven different reagent ions generated by the electrospray. We analyzed the mobility-mass spectra of the oxidation products C9-10H14-18O2-6. Our results show that acetate and chloride yield the highest charging efficiencies. Analysis of the mobility spectra suggests that the clusters have 1-5 isomeric structures (i.e., ion-molecule cluster structures with distinct mobilities), and the number is affected by the reagent ion. Most of the isomers are likely cluster isomers originating from binding of the reagent ion to different sites of the molecule. By comparing the number of observed isomers and measured mobilities and collision cross sections between standard pinanediol and pinonic acid to the values observed for C10H18O2 and C10H16O3 produced from oxidation of α-pinene, we confirm that pinanediol and pinonic acid are the only isomers for these elemental compositions in our experimental conditions. Our study shows that the SESI-DMA-MS produces new information from the first steps of oxidation of α-pinene.

Gas-Phase Peroxyl Radical Recombination Reactions: A Computational Study of Formation and Decomposition of Tetroxides.

The recombination ("dimerization") of peroxyl radicals (RO2•) is one of the pathways suggested in the literature for the formation of peroxides (ROOR', often referred to as dimers or accretion products in the literature) in the atmosphere. It is generally accepted that these dimers play a major role in the first steps of the formation of submicron aerosol particles. However, the precise reaction pathways and energetics of RO2• + R'O2• reactions are still unknown. In this work, we have studied the formation of tetroxide intermediates (RO4R'): their formation from two peroxyl radicals and their decomposition to triplet molecular oxygen (3O2) and a triplet pair of alkoxyl radicals (RO•). We demonstrate this mechanism for several atmospherically relevant primary and secondary peroxyl radicals. The potential energy surface corresponds to an overall singlet state. The subsequent reaction channels of the alkoxyl radicals include, but are not limited to, their dimerization into ROOR'. Our work considers the multiconfigurational character of the tetroxides and the intermediate phases of the reaction, leading to reliable mechanistic insights for the formation and decomposition of the tetroxides. Despite substantial uncertainties in the computed energetics, our results demonstrate that the barrier heights along the reaction path are invariably small for these systems. This suggests that the reaction mechanism, previously validated at a multireference level only for methyl peroxyl radicals, is a plausible pathway for the formation of aerosol-relevant larger peroxides in the atmosphere.

Energy transfer, pre-reactive complex formation and recombination reactions during the collision of peroxy radicals.

In this paper we study collisions between polyatomic radicals - an important process in fields ranging from biology to combustion. Energy transfer, formation of intermediate complexes and recombination reactions are treated, with applications to peroxy radicals in atmospheric chemistry. Multi-reference perturbation theory, supplemented by coupled-cluster calculations, describes the potential energy surfaces with high accuracy, including the interaction of singlet and triplet spin states during radical recombination. Our multi-reference molecular dynamics (MD) trajectories on methyl peroxy radicals confirm the reaction mechanism postulated in earlier studies. Specifically, they show that if suitable pre-reactive complexes are formed, they will rapidly lead to the formation and subsequent decomposition of tetroxide intermediates. However, generating multi-reference MD trajectories is exceedingly computationally demanding, and we cannot adequately sample the whole conformational space. To answer this challenge, we promote the use of a novel simplified semi-empirical MD methodology. It assumes the collision is governed by two states, a singlet (S0) and a triplet (T1) state. The method predicts differences between collisions on S0 and T1 surfaces, and qualitatively includes not only pre-reactive complex formation, but also recombination processes such as tetroxide formation. Finally, classical MD simulations using force-fields for non-reactive collisions are employed to generate thousands of collision trajectories, to verify that the semi-empirical method is sampling collisions adequately, and to carry out preliminary investigations of larger systems. For systems with low activation energies, the experimental rate coefficient is surprisingly well reproduced by simply multiplying the gas-kinetic collision rate by the simulated probability for long-lived complex formation.

Odd-Number Cyclo[n]Carbons Sustaining Alternating Aromaticity.

Cyclo[n]carbons (n = 5, 7, 9, ..., 29) composed from an odd number of carbon atoms are studied computationally at density functional theory (DFT) and ab initio complete active space self-consistent field (CASSCF) levels of theory to get insight into their electronic structure and aromaticity. DFT calculations predict a strongly delocalized carbene structure of the cyclo[n]carbons and an aromatic character for all of them. In contrast, calculations at the CASSCF level yield geometrically bent and electronically localized carbene structures leading to an alternating double aromaticity of the odd-number cyclo[n]carbons. CASSCF calculations yield a singlet electronic ground state for the studied cyclo[n]carbons except for C25, whereas at the DFT level the energy difference between the lowest singlet and triplet states depends on the employed functional. The BHandHLYP functional predicts a triplet ground state of the larger odd-number cyclo[n]carbons starting from n = 13. Current-density calculations at the BHandHLYP level using the CASSCF-optimized molecular structures show that there is a through-space delocalization in the cyclo[n]carbons. The current density avoids the carbene carbon atom, leading to an alternating double aromaticity of the odd-number cyclo[n]carbons satisfying the antiaromatic [4k+1] and aromatic [4k+3] rules. C11, C15, and C19 are aromatic and can be prioritized in future synthesis. We predict a bond-shift phenomenon for the triplet state of the cyclo[n]carbons leading to resonance structures that have different reactivity toward dimerization.

Pathways to Highly Oxidized Products in the Δ3-Carene + OH System.

Oxidation of the monoterpene Δ3-carene (C10H16) is a potentially important and understudied source of atmospheric secondary organic aerosol (SOA). We present chamber-based measurements of speciated gas and particle phases during photochemical oxidation of Δ3-carene. We find evidence of highly oxidized organic molecules (HOMs) in the gas phase and relatively low-volatility SOA dominated by C7-C10 species. We then use computational methods to develop the first stages of a Δ3-carene photochemical oxidation mechanism and explain some of our measured compositions. We find that alkoxy bond scission of the cyclohexyl ring likely leads to efficient HOM formation, in line with previous studies. We also find a surprising role for the abstraction of primary hydrogens from methyl groups, which has been calculated to be rapid in the α-pinene system, and suggest more research is required to determine if this is more general to other systems and a feature of autoxidation. This work develops a more comprehensive view of Δ3-carene photochemical oxidation products via measurements and lays out a suggested mechanism of oxidation via computationally derived rate coefficients.

Magnetically induced ring currents in metallocenothiaporphyrins.

The magnetically induced current-density susceptibility tensor (CDT) of the lowest singlet and triplet states of the metallocenothiaporphyrins, where the metal is V, Cr, Mn, Fe, Co, Ni, Mo, Tc, Ru, or Rh, have been studied with the gauge-including magnetically induced currents (GIMIC) method. The compounds containing V, Mn, Co, Tc or Rh were studied as cations because the neutral molecules have an odd number of electrons. The calculations show that the aromatic nature of most of the studied molecules follows the Hückel and Baird rules of aromaticity. CDT calculations on the high-spin states of the neutral metallocenothiaporphyrins with V, Mn, Co, Tc or Rh also shows that these molecules follow a unified extended Hückel and Baird aromaticity orbital-count rule stating that molecules with an odd number of occupied conjugated valence orbitals are aromatic, whereas molecules with an even number of occupied conjugated orbitals are antiaromatic.